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Abstract 

An impurity bond J' in a periodic ID antiferromagnetic, spin 1 chain with 

exchange J is considered. Using the numerical density matrix renormahza- 

tion group method, we find an impurity energy level in the Haldane gap, 

corresponding to a bound state near the impurity bond. When J' < J the 

level changes gradually from the edge of the Haldane gap to the ground state 

energy as the deviation dev = ( J — J')/ J changes from to 1. It seems that 

there is no threshold. Yet, there is a threshold when J' > J. The impurity 

level appears only when the deviation dev = (J' — J)/ J' is greater than 3^ 

which is near 0.3 in our calculation. 
PACS numbers:75.10Jm,76.50+g,75.50.Ee 
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I. INTRODUCTION 



The Heisenberg model of ID antiferromagnetic(AF) chains is : 



H = jj:Si-s 



•i+l 



(1) 



with J > 0, and Sj is a quantum spin. Up till now, the exact solution of this Hamiltonian 
is available only when s = 1/2. The solution is obtained by the Bethe ansatz which 
shows there is no gap in the energy spectrum. Then, what is the situation when s = 1 or 
other integer? Some years ago, Haldane made his famous conjecture that ID AF chains 
of integer spin have a gap, while those of half integer spin have a gapless spectrum. Since 
then, a lot of work has been done. In particular, a Valence-Bond-Solid (VBS) picture was 
proposed by Affleck, Kennedy, Lieb, and Tasaki(AKLT) [|,|| to interprete the ground state 
of the integer spin AF chains. The proposed picture agrees quite well with both experimental 
P and numerical studies. 

An important issue of integer spin AF chains is the doping effect which involves some fun- 
damental many-body quantum problems . DiTusa et al p| have done doping experiments 
with Zn or Ca in Y2BaNiO^, which contains ID AF chains of Ni — (O) ~ Ni — [O) — ■ ■ -. 
In the Zn doping case, the non-magnetic Zn^^ ion substitutes Ni^^ to sever the AF chain 
giving rise to finite length effect. On the other hand, in the Ca doping case {Ca^^ substitut- 
ing the off-chain atom Y^^), holes are introduced in the oxygen orbitals along the Ni — O 
chain, modifying the superexchange interaction and producing an impurity state inside the 
Haldane gap, as seen by neutron scattering experiments. The latter case can be represented 
by changing J to J' at the impurity bond. We focus on this case in this Brief Report. Thus 
we write down our Hamiltonian as: 



where the impurity bond is put between sites and 1. 

Recently a calculation of the dynamic structure factor based on the Schwinger-boson 
approach was performed for the bond-doping case which first indicated the existence of 
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if = J X) Sj ■ Sj+i + (J' — J)So ■ Si, 



(2) 



a threshold for the bond-couphng deviation. According to this calculation, a bound state 
is induced by the impurity only if the deviation exceeds this threshold, and the impurity 
level is located in the middle of the Haldane gap, almost independent of bond-coupling when 
temperature is much lower than the Haldane gap. 

The recently developed numerical method of density matrix renormalization group 
(DMRG) P,|lOi has achieved a great success in calculating the low energy spectrum of ID 
Heisenberg AF chains so it's interesting and natural to use this method to study the 



doping effect on the Haldane gap systems. Such a calculation was carried out by Sorenson 
and Affleck (SA) They reported the existence of the threshold, but the impurity en- 
ergy level they obtained changes gradually as the strength of impurity bond changes. The 
existence of the threshold is in agreement with the Schwinger boson calculation [0], while 
the gradual change of the impurity level is in disagreemnet with the latter. However, these 
authors have used the open boundary conditions and the way they added new sites may 
affect the bound state itself (See the discussion in the next Section). To clarify this issue we 
have decided to devise a more elaborate scheme to carry out the DMRG calculations. 

II. DMRG METHOD 



The DMRG method has been developed by S.R. White et al P,p!0| in recent years. 
Generally, the standard pattern looks like this: 



* • ■ ■ * — 


— * — 


— * — 





During each iteration, the chain (superblock) is divided into 4 blocks, where the left- 
most and the right-most (which can be named as old blocks) are the same for the symmetry 
consideration. The two middle blocks contain only two sites. After each run, the left and 
right halves of the superblock are optimized to form new blocks which are used for the next 
iteration. This is, however, not good in the presence of an impurity bond. At first glance 
it seems right to put the impurity bond in the middle of the chain (this is what SA have 



done |T2|), but this configuration may introduce some artificial effects. For example, the 
two middle sites which were optimized with the impurity bond between them, are used in 
later calculations without impurity being close to them. When the chain is very long, the 
impurity bond effect should go to zero during iterations (since there is only one impurity 
bond), but using this pattern, the effect cannot become small as iterations go on. 

To calculate the impurity level of the periodic chain, we do our numerical work using the 
following pattern: 




The impurity bond is placed at the left-hand side, and we add two sites each time with 
one each at or near the center of the upper and lower parts of the above periodic chain. 
Using this pattern, one can avoid the difficulties mentioned above. We preserve m = 100 
states in each run and discard the rest. Since Sz is a good quantum number, we use it to 
reduce our Hamiltonian matrix. 

Shank transformation has been demonstrated as a powerful extrapolation method to 
get Haldane gap precisely |T^. We use it to extrapolate our results to the infinite length 



limit. To verify the credibility of this method, we also use Aitkens' transformation [0 to 
extrapolate our finite chain data. They both give very similar values (usually the difference 
is less than 10^"^). Compared with the previous conclusion in special cases, our results show 
reasonable precision. So we believe that our value is correct at least up to 10~^. 



III. NUMERICAL RESULT 



We use VBS picture to give some intuitive explanation of our calculation results, 
although this picture is not quite precise. According to this picture, when considering the 
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low excitation energy spectrum, the spin 1 of each site can be regarded as two spin 1/2's 
combined together to form a triplet S = 1 state, yet the two 1/2 spins of the nearest sites 
form a Valence-Bond, or, a singlet state. Therefore, the ground state of an open chain can 
be labelled by the sum of the two edge sites 1/2 spins. The ground state of the ID AF chain 
is four-fold degenerate, with one state of Stotai = 0, where the two 1/2 spins of each edge 
form a singlet, and three states of Stotai = 1 , where the two 1/2 spins form a triplet. The 
first excited states are the quintlet states with Stotai = 2, which lie above the ground state 
by the famous Haldane gap. Above them, the continuous spectrum starts. For a periodic 
chain, VBS picture also holds, but the triplet state rises up to the edge of the continuous 
energy spectrum, leaving the ground state nondegenerate. 

The impurity state of the periodic chain, is just the Stotai = 1 triplet state. When 
J'/ J = 0, the periodic chain becomes an open chain, and this triplet will merge with the 
Stotai = ground state in the infinite chain length limit. Thus the impurity energy level is 
at 0. On the contrary, when J'/ J = 1, the periodic chain recovers, so the "impurity" level 
will rise up to merge with the Haldane gap. 

Because of the SO (3) symmetry of the Heisenberg Hamiltonian , the Stotai = 1 triplets 
are degenerate, so we choose the state = 1 as a representative of them. Also we set the 
zero of the energy spectrum as the ground state energy. 

First we consider the case < J' < J, which means the super- exchange is smaller due to 
doping, yet it does not change the AF nature. The calculated result is presented in Figure 1, 
from which we see that the impurity level increases linearly near dev = l{dev = {J — J') / J) , 
which is in accordance with SA's perturbation argument |T^. However, compared with SA's 



work |T2[ , our numerical result favors the conclusion that there is no threshold. As seen from 
the figure, the impurity energy appears as soon as dev ^ 0, and shows a quadratic behavior 
near dev = 0. Kaburagi et al [|T5| has investigated the impurity bond effect in terms of 
domain wall exitations. Their variational results show the linear behavior near dev = and 
quadratic behavior near dev = 1. They also perform diagonalizations by Lanczos method 
for = 12 and 13 chains to support the variational results. 



When J' > J, the result is shown in Figure 2, where the impurity level vs dev = 
{J' — J)/ J' is plotted. This time, a threshold B^. appears near 0.3, i.e., only when dev > Be 
an impurity localized state will appear. Again the level rises linearly near dev = 1, in 
agreement with the results of SA [|12|, but with a different Be (they gave Be = 0.5). The 
case when J'/ J oo {dev 1) can be understood easily as follows: The two spins of the 
impurity bond form a singlet state because of the strong coupling, breaking the periodic 
chain of length L to make an open chain of length L-2. So, as following from our calculation, 
the difference between the impurity level and the ground state will go to zero in this limit. 
According to the finit chain results, the impurity level when < dev < Be is greater than 
the level when dev = 0, and the difference between them become smaller when the chain 
length increases, so we believe there is a threshold when the length goes to infinity, instead 
of the exponential or other behavior of the energy near dev = . 

Moreover, to ensure that the impurity level means a local state near the impurity bond, 
we plot < Szi > of the periodic chain in Figure 3 (J' < J). When J'/ J = 0{dev = 1), the 
"boundary 1/2 spin" of open chain appears, just as the VBS picture assumes [Q]. When dev 
goes from 1 to 0, the state becomes more and more delocalized. Figure 4 shows the case of 
J' > J, where the two sites at the middle form a singlet state, while the sites near them 
show the "boundary 1/2 spin", too, when dev 1. 

Figure 5 and 6 give the Shank transformation results of \dE/dDev\. There is obviously 
no abrupt change. 

IV. CONCLUSION 

Using an elaborate numerical DMRG method, we have investigated the bond doping 
effect of ID AF chain. A threshold Be exists only when J' > J. When J' < J ot J' > J 
but dev < Be, no impurity state is induced. When dev > Be and J' > J, the impurity 
level within the Haldane gap corresponds to a local state near the impurity bond, and the 
level changes gradually when dev runs from Be to 1. The discrepancy of the present study 
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with the Schwinger boson calculation [0 is not fully understood at present. One possiblity 
is that the Schwinger boson representation introduces an additional symmetry (which sets 
the impurity level at the middle of the Haldane gap), not inherent to the original quantum 
spin systems. 

Our calculation is supported by LSEC, COAST and Laboratory of Oomputational 
Physics, lAPCM. We appreciate their assistance. After completion of our calculations we 
received a preprint by Xiaoqun Wang and S. Mallwitz, which reports essentially the same 
results. We thank Dr. Wang for sending us the preprint. 
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FIGURES 

FIG. 1. Impurity energy level vs impurity dev = ( J — J')/ J, when < J' < J. The chain is 
periodic, and the dashed line indicates the Haldane gap. 

FIG. 2. Impurity energy level vs impurity dev = (J' — J)/ J', when J' > J. The chain is 
periodic, and the dashed line indicates the Haldane gap. 

FIG. 3. < Szi > of the impurity state as function of chain index i when J' < J. The 
star is for the case of dev = 1, while the triangle is for dev = 0.6, and the square for dev = 
{dev = ( J — J')/ J). The impurity bond is placed between site 15 and 16. 

FIG. 4. < Szi > of the impurity state as funtion of chain index i when J' > J. The star 
is for the case of dev = 0.999, while the triangle is for dev = 0.7, and the square for dev = 
{dev = {J' — J)/ J'). The impurity bond is placed between site 15 and 16. 
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